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ABSTRACT 



Context. The nature of Type la supernova progenitors is still unclear. The outstanding characteristic of the single-degenerate scenario 

is that it contains hydrogen in the binary companion of the exploding white dwarf star, which, if mixed into the ejecta of the supernova 

in large amounts may lead to conflicts with the observations thus ruling out the scenario. 

Aims. We investigate the effect of the impact of Type la supernova ejecta on a main sequence companion star of the progenitor system. 

With a series of simulations we investigate how different parameters of this system affect the amount of hydrogen stripped from the 

companion by the impact. 

Methods. The stellar evolution code GARSTEC is used to set up the structure of the companion stars mimicking the effect of a binary 

evolution phase. The impact itself is simulated with the smoothed particle hydrodynamics code GADGET2. 

Results. We reproduce and confirm the results of earlier grid-based hydrodynamical simulation. Parameter studies of the progenitor 

system are extended to include the results of recent binary evolution studies. The more compact structure of the companion star found 

here significantly reduces the stripped hydrogen mass. 

Conclusions. The low hydrogen masses resulting from a more realistic companion structure are consistent with current observational 

constraints. Therefore, the single-degenerate scenario remains a valid possibility for Type la supernova progenitors. These new results 

are not a numerical effect, but the outcome of different initial conditions. 

Key words, stars: supernovae: general - hydrodynamics - binaries: close 



1. Introduction 

While the progenitors for Type Ib/c and Type II super- 
novae are known. Type la supernovae (SNe la) still elude 
an identification of their progenitor system. This is an un- 
pleasant situation given the fact that these objects are one 
of the most important tools to determine cosmological pa- 
rameters. By virtue of empirical calibration methods (e.g.. 



iPhillipsI Il993h they can be used as standardizable candles 
for distance measurements. This calls for an understanding 
of the mechanism of SNe la; and indeed, some progress 
has been made in recent years in understanding the explo- 
sion mechanism in terms of thermonuclear expl o sions of white 



'Reinecke et al ' 



dwarf (WD) stars (e.g 

200l iRopke & Hi llebrandt 2005 



2002 : Gamezo et al.l 

, Rop ke & Niemeveij 120071 

Mazzah et al.ll2007l; iRopke et al.ll2007h . In order to judge poten- 



tial systematic errors in SN la cosmology, a theoretical connec- 
tion between the explosion characteristics and properties of the 
progenitor system would be desirable. Yet despite all efforts on 
both the theoretical and on the observational side, the nature of 
the progenitor system remains enigmatic. 

The stabilization of WDs against gravity does not depend on 
a finite energy source such as the nuclear burning in normal stars. 
Due to the Fermi pressure of a degenerate electron gas, single 
WD stars are in principle eternally stable. Thus, some additional 
dynamics is required in order to reach an explosive state. The 
most likely possibility is a WD being part of a binary system 
and accreting matter from its companion. 

Current progenitor models distinguish between the single 
degenerate and the double degenerate scenario. The former 



(proposed bv IWhelan & Ibeiilll973l) assumes a "normal", non- 
degenerate star to be the binary companion - either a main se- 
quence (MS) star or a red giant (RG). In this case, the WD 
accretes mass from its MS or RG companion via Roche-lobe 
overflow or by winds (symbiotic systems) until it approaches 
the Chandrasekhar mass. The densities reached in the core of 
the WD are then sufficiently high to trigger nuclear reactions 
which finally cause a thermonuclear explosion of the star (but 
note that an explosion be fore reaching the Chandrasekhar mass 
may also be possible, e.g. Fink et all2007l). The do uble degener- 
ate scenario (llben & Tutukovll 984 IWebbinkI 19841) . on the other 
hand, assumes a binary system of two WDs. Due to gravitational 
wave emission, the system becomes unstable at some point, the 
WDs merge, and may eventually explode in a SN la. A sum- 
mary of arg ument s in fav or of and against both scenarios can 
be found in iLiviol (|2000|) . In theoretical modeling, the single- 
degenerate Chandrasekhar-mass scenario has received most at- 
tention recently. By constraining the amount of fuel available in 
the thermonuclear explosion, it provides a natural explanation 
for the observed uniformity of SNe la. 

Observationally, its hard to distinguish between the progen- 
itor scenarios. One fundamental difference is the complete ab- 
sence of hydrogen in the double degenerate scenario as it as- 
sumes a merger of two carbon/oxygen WDs. In contrast, in the 
typical single degenerate scenario, hydrogen is the main con- 
stituent of the companion. An exception are helium-accretors in 
which a more evolved companion star has lost his hydrogen en- 
velope. Thus, the WD accretes heli um instead of hydrogen and 
hydrogen is missing in the system. iKato & Hachisul (l2003l) re- 
ported a possible detection of such an object. In the standard 
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scenario, however, the companion star features a hydrogen enve- 
lope; and at least some part of it is expected to be carried away 
by the SN la ejecta impacting the companion. This, in princi- 
ple, causes a problem for the single-degenerate scenario, because 
the astronomical classification of SNe la rests on the absence of 
hydrogen-features in the spectra of these events. The hydrogen 
stripped off from the companion will have rather low velocities. 
It may thus be detectable in nebular spectra, if abundant enough. 
For the single-degenerate scenario it is therefore of critical im- 
portance that the mass of stripped material is sufficiently low to 
be still consistent with the observations. 

There have been a few att empts to search for hydrogen in 
nebular spectra of SN la. Mattila et al.l ( 120051) studied nebular 
spectra of SN 200 lei. From modeling them, they derived an up- 
per limit of 0.03 M© of solar abu ndance material at velocities 
lower than lOOOkms '. Recentlv. [Leonard (l2007h studied neb- 
ular spec tra of SN 2005am a nd SN 2005cf. Based on the same 
model as lMattila et al.l (l2005b . he estimated <0.01 Mq of hydro- 
gen material for both objects. Recently also hydrogen has been 
detected indirectly by iPatat et aTl ( |2007[) in circumstellar mate- 
rial of SN 2006X. 

An alternative to this approach of observationally constrain- 
ing the nature of the progenitor system is to directly search for 
the former companion star of the single-degenerate scenario in 
the remnants of historical galactic SNe la. Such a search has 
been carrie d out in the remnan t of Ty cho Brahe's supernova 
of 1572 bv iRuiz-Lapuente et al.l (|2004|) . who claimed the iden- 
tification of the binary companion. The star in question is a 
slightly evolved solar-type star, that moves with a radial velocity 
of -108 km s"' relative to the sun. It also has an atypical large 
tangential velocity of about 90 km s"' . A significantly larger ve- 
locity of the star compared to neighbours is expected as a result 
of the disappearing binary orbit. Other stars observed in the same 
area with similar distances move only with average radial veloc- 
ities of about -20 to -40 km s"', with a velocity dispersion of 
about 20km s"'. 

On the theory side, [Marietta et al.l ( l2000l) presented two- 
dimensional hydrodynamical simulations of the impact of 
SNe la on their companions. They found that 0.15 M© were 
stripped from a Roche-lobe filling MS companion. Th is would 
rule o ut a MS-WD system for the SN la a nalysed bv iLeonardI 
(l2007l) . if it was representative. Recendy, iMeng et al.l (l2007l) 
pointed out that considering the effect of the mass transfer phase 
on the companion star may change the result significantly. They 
studied the impact of SN e la on diff'ere n t com panion stars an- 
alytically. In contrast to [Marietta et al.l (|2000|) . w ho assumed 
the str ucture of single MS stars for the companion, iMeng et all 
(|2007|) evolved it through the binary evolution phase before the 
explosion of the WD. They found at least 0.035 M© of stripped 
hydrogen for the companion. However, this result is only a 
lower limit, since they did not include mass loss by vaporiza- 
tion from the hot surface of the star. Thus, taken at face value, 
the currently available theoretical studies constitute a strong case 
against MS-hWD progenitor systems for SNe la. 

The a im of the study prese nted here is to check and up- 
date the [Marietta et al.l (l2000h calculations with the results 
of rec ent detailed binary evolution models. Ilvanova & TaamI 
(|2004|) identified possible SN la progenitors from a parameter 
study of MS+WD binary evolution. There results are in agree- 
ment with other studies in this field (e.g. JLanger et al.l l2000t 
iHan & Podsiadlowskill2004b . Based on these results, we present 
an exploration of the effect of the impact of a SN la on differ- 
ent MS companions by 3D hydrodynamical simulations. Section 
|2] summarizes the codes used. Section |3] demonstrates that our 



approach reproduces the results of [Marietta et al.l ( l2000l) and 
presents a resolution study. Section|4]discusses an exploration of 
different progenitor systems and Section |5]derives observational 
implications of our results. A summary and an outlook conclude 
this work in Section|6] 

2. Modeling approach 

Two different codes are employed in this work: one to construct 
the companion stars mimicking a binary evolution and the other 
to investigate the hydrodynamical impact of the supernova on 
the companion star. 

To evolve the co mpanion stars we use th e stellar evolution 
code GARSTEC of IWeiss & Schlattll (|2007|) . It evolves stars 
with a given mass and metallicity to a certain age and is used to 
constru ct a solar-type compa nion star similar to the "HCV" sce- 
nario of [Marietta et al.l (|2000|) (see Sect.|3]l. The code also allows 
to include mass loss during the evolution but does not account 
for a binary evolution. In our study of a variety of progenitor 
models (see Sect. |4]i, we t herefore rely on the para meters of the 
binary evolution study by Ilvanova & TaamI ( 120041) to construct 
our companion stars. For each of the models we first set up a star 
that fits the parameters at the onset of the mass transfer phase. 
At this point, a constant mass loss rate is assumed and the stellar 
evolution is followed for the duration of the mass transfer period. 
The mass loss rate employed here corresponds to the mass loss 
rate of the original binary models of Ilvanova & TaarnI ( |2004|) av- 
eraged over the entire mass transfer phase. After following this 
phase, we obtain a stellar configuration that approximates the 
outcome of a realistic binary evolution. It is used to study the 
process of the interaction with the explosion ejecta once the bi- 
nary WD undergoes a SN la. 

The impact of the supernova ejecta is simulated using the 
smoothed parti cle hydrodynamics (SPH) code GADGET2 of 
ISpringell (l2005h . In order to set up the companion star here, we 
map the one-dimensional profiles of density, internal energy, and 
nuclear composition of the stellar evolution calculation to a par- 
ticle distribution suitable for the SPH code. The mapping is done 
by transforming a uniform particle distribution to the given ra- 
dial density profile. Details of this procedure wiU be explaind in 
(iPakmor et al.ll2008l) . Using GADGET2, these companion stars 
were relaxed in a separate step for 1 .0 x 1 0"* s to get rid of numeri- 
cal artifacts that may have been introduced by the transformation 
(e.g. due to the random placement of the particles). 

A supernova was added to the simulation at a distance given 
by the orbital period of the binary system before the explo- 
sion. The supernova was set up based on the W7 model by 
iNomoto et al.l(ll984h . This one-dimensional model is well tested 
and provides a good fit to observations of standard SN la. It has 
a kinetic energy of 1.23 x 10^' erg. At the time we add the super- 
nova model to the simulation, it has reached already the phase 
of homologous expansion. The impact of the supernova ejecta 
and the following evolution of the ejecta and the companion star 
are simulated for about one hour After this time, the companion 
star is already relaxing and its mass and velocity have reached 
constant values. 

The GADGET code was so far on ly used for cosmological 
simulations , with few exceptions (e.g.. lMorris & Podsiadlowskil 
l2006il2007l) . To be able to use it for a stellar astrophysics prob- 
lem, a few extensions of the original implementation were nec- 
essary. The following basic setup was used in these simulations: 

- The smoothing length is chosen such that a sphere of its ra- 
dius encloses 80 neighboring particles. 
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- The gravitational softening length equals the smoothing 
length. 

- All particles are given the same mass. 

- Nuclear reacti ons are neglec t ed in the simulations. This is 
justified since [Marietta et al.l (|2000|) showed that the addi- 
tional energy generated by shock-wave induced hydrogen 
burning is a marginal efl'ect. The savings in computing time 
due to this approximation, however, are substantial. 

Details of these cha nges are described in a separate publication 
(lPakmoretal.1120081) . 



3. Tests of implementation 

One of the obvious questions arising in our approach is whether 
the 3D SPH scheme applied here leads to the same results as 
the 2D grid-based approach of Marietta et al.] (|2000|) . This is 
tested by using the initial parameters of the HCV scenario of 
[Marietta et al.l ( 120001) in our setup. In the HCV scenario, the su- 
pernova is reahzed as a W7 model. The companion is a solar like 
1 .017 Mo main sequence star with a central hydrogen abundance 
of 0.58. The separation between supernova and companion star 
at the time of the explosion is 2.04 x 10" cm. 

Figure [T] shows the typical evolution of the companion star 
in our simulations. Here, an example with a total of 235499 SPH 
particles is illustrated starting out with the impact of the outer- 
most supernova ejecta on the companion and ending when they 
have passed the star and it relaxes again. The first image af- 
ter 20 s shows the companion star at the instant when the first 
ejecta reach the companion from the right. In the second snap- 
shot, taken after 200 s, the ejecta have hit the companion star. 
A shock wave forms and starts to propagate through it. Another 
200 s later (third snapshot of Fig. [U, the shock wave reaches the 
center of the star In the fourth snapshot (~2000 s after the explo- 
sion the shock wave has crossed the star completely. Material is 
ejected on its far side. The last two images show the star shrink- 
ing and relaxing again. 

Qualitatively, th is looks similar to the simulations of 
[Marietta et al.l ( l2000l) . However, our results do not show the hy- 
drodynamical instabilities they observed at the interface of su- 
pernova ejecta and the stripped material in the wake of the star 
(note that some mixing of companion star material into the su- 
pernova ejecta is realized, see Fig. [TJ. This difference is not sur- 
prising as SPH codes are known to suppress i nstabilities due to 
their large numerical viscosity (see, however. iFrver et al1l2007l 
showing that grid-based codes may under some circumstances 
not reproduce the mixing better than SPH codes). However, it is 
not a priori clear that this morphological difference significantly 
affects the effects and quantities we are interested in. The fun- 
damental quantity upon which we base our comparison is the 
bound mass of the companion star, or, equivalently, the mass 
stripped away from it by the supernova ejecta. 

As a first step, we p erform a resolu t ion stu dy in order ensure 
a comparison with the [Marietta et al.l ( [20001) result based on a 
converged simulation. To test for numerical convergence we use 
different resolutions with 1.1 x 10^ to 4.7 x 10^ SPH particles 
(counted for the entire setup, i.e. supernova and companion star, 
as the total mass of the configuration is equally distributed on 
the particles). 

Figure|2]shows the evolution of the mass loss from the com- 
panion star with time for different resolutions. This value is cal- 
culated by subtracting the sum over the masses of all particles 
that are gravitationally bound to the star from its initial mass. 
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Fig. 2. Mass loss of the companion star depending on the time 
after the supernova explosion for different resolutions. 

Whether or not a particle is bound to the star is decided by com- 
paring its potential energy relative to the center of the star with 
the kinetic energy of its motion relative to the motion of the star 
For this, the center and velocity of the star are taken from the 
previous snapshot. Then the current center and velocity of the 
star are calculated from all bound particles. In principle, the new 
position of the star's center should be used to recalculate which 
particles are bound to the star and this cycle should be iterated 
until convergence is reached. However, already the position and 
velocity of the star taken from the previous snapshot usually pro- 
vide a sufficiently accurate approximation, and we thus forego 
the iteration. 

After ~3000s the mass loss has stopped and the star has 
reached its final mass. There is a numerical artefact in the de- 
tection of the bound mass of the companion star around ~300 s 
after the explosion. For one snapshot the mass of the companion 
star seems to increase. This is due to a false detection of particles 
as unbound in the previous snapshot due to high radial velocities 
of particles, that point to the center of the star 



"star 


«tot 


'"particle 


'^Stripped 


'^^ remnant 


"kick 






[Mo] 


[Mo] 


[Mo] 


[kms-'] 


50000 


117738 


2.03 x 10-5 


0.126 


0.891 


101.1 


100000 


235499 


1.02 x 10-5 


0.144 


0.873 


95.4 


500000 


1177482 


2.03 X 10-* 


0.138 


0.879 


85.3 


1000000 


2352965 


1.02 X 10-* 


0.135 


0.882 


80.3 


2000000 


4709926 


5.09 X 10-' 


0.134 


0.884 


81.6 



Table 1. HCV resolution test: «star, «tot, and ntparticie denote the 
number of particles the star is composed of, the total number of 
particles in the simulation, and the particle mass, respectively. 
The results of the simulations are characterized by mass Wstripped 
stripped away from the companion star, its final mass m,-emnant, 
and its velocity relative to the supernova i>kick- All values are 
taken 2 x 10"* s after the explosion. 



Table [T] shows again some simulation properties for differ- 
ent resolutions after 2 x 10"* s. Considering Table [1] and Fig.|2] 
we conclude, that this simulation is numerically converged when 
using more than 10^ particles. For the two simulations with 2.3 
and 4.7 million particles, the graphs in Fig.|2]are nearly identical. 
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Fig. 1. Snapshots of the evolution of the companion star in the HCV scenario. The plots use cylindrical coordinates. The radial 
coordinate is averaged over angle. Color-coded is the density. 



Note that the kick velocity converges slowlier than the stripped O. 134M0. This is quite close to the result of 0.15 M© reported 

mass. bv [Marietta et al.1 (l2000h . 

As carried out here, however, the c omparison is n ot ye t 

Numerical convergence, however, does not necessarily im- based on exactly the same assumptions. [Marietta et al.l (|2000|) 

ply consistency with the physical solution. To check our re- performed their simulations on a finite computational grid and 

suits in this r e spect, we compare them to previous results of therefore mass was lost over the domain boundaries. This mass 

[Marietta et al.l ( l2000h . With the highest resolution of nearly was always assumed to be unbound from the companion star 

5 million particles the stripped mass after 2.0 x IC^s is We therefore recalculate the stripped mass in our simulations 
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by assuming for comparison all particles to be unbound that 
are outside a cylindrical box equal to the simulation volume of 
[Marietta et al. (2000). In this approach we find a stripped mass 
of 0.143 M0 for our highest resolved simulation. This result is 
in excellent agreement w ith the 0.15 M© of stripped mass found 
bv [Marietta et al.l (|2000|) . Moreover , the remnant star ve locities 
at this time (85.7 km s"' reported bv lMarietta et alj|2000l vs. our 
81.4kms"') agree very well. 

We therefore conclude that our SPH approach is capable of 
capturing the main dynamical effects of the supernova i mpact on 
the co mpanion star. The global quantities of the Marie tta et alj 
(I2OOOI) study are reproduced down to the percent level. The dif- 
ferences in the occurrence of hydrodynamic instabilities are ob- 
viously a minor effect with respect to the overall results such as 
the stripped mass and the velocity of the companion star caused 
by the kick by the supernova ejecta. The instabilities suppressed 
in our SPH approach may, however, enhance the mixing between 
supernova ejecta and stripped material from the companion in 
reality. 



4. 1 . Explosion energy 

The influence of the supernova explosion energy on the inter- 
action with the companion is studied on the basis of model 
rp3_20a. All parameters but the supernova energy were kept con- 
stant with the values of the original model (see Table|2li. 

The kinetic energy of the supernova was varied in the 
range 0.8 ... 1.6 B (= 0.8 ... 1 .6 x 10^' erg). The lower limit cor- 
responds to lowest kinetic energies of simulated deflagration 
SNe la. The upper limit is the maximum energy a SN la can 
have, assuming that a Chandrasekhar-mass WD consisting of an 
equal-by-mass mixture of C and O burns completely to Ni^*. 
The kinetic energy of the supernova ejecta E^^^j^ was adjusted 
by scaling the velocities v' of the supernova particles (originally 
representing the W7 model with E*^„j^ and v^'') according to 



A 



F' 

kin.SN W7 
■ ■ V . 

pW7 

■"kin.SN 



(1) 



4. Parameter studies 



This scaling preserves the homologous expansion (u oc r) of the 
ejecta. 



Three major physical parameters of the progenitor system are 
expected to influence the dynamics of the supernova impact on 
the companion star: the kinetic energy of the supernova ejecta 
(powered by the thermonuclear burning in the explosion), the 
separation between supernova and companion, and structure of 
the companion star at the time of the explosion. 

In reality, of course, separation and companion structure are 
not independent, but result from the characteristics of the orig- 
inal binary system and its evolution through the mass transfer 
phase. Moreover, the dependence of the supernova explosion en- 
ergy on other parameters is unknown and thus treated as an in- 
dependent physical parameter in this study. In any system with 
Roch-lobe overflow, the solid angle under which the companion 
star is seen from the white dwarf depends only on the mass ratio 
of its components. However, in order to determine how changing 
the distance affects the results it is treated here as an independent 
parameter. This is motivated partially by the possibility that the 
WD and the compani on star disconnect shortly before the explo- 
sion, as suggested bv lPatat et all (l2007h for SN 2006X. 

Our main emphasis, however, will be on the structure of 
the companion star. In order to model the progenitor systems , 
we take parameters from the study of llvanova & TaamI (|2004|) 
and a typical Populat ion I metallicity of 0.02 for all stars. 
llvanova & TaamI ( |2004|) analyzed the evolution of possible SN la 
progenitor systems consisting of a WD and an evolved MS star 
Detailed calculations of 65 different systems were carried out, 
varying the initial mass of both objects and the initial distance 
to cover the parameter space. From the resulting sample, they 
pick six representative models which are believed to be likely 
to evolve into a SN la. These are the systems we select for our 
study. Their properties are listed in table |2| 

The three parameters are individually discussed in the fol- 
lowing paragraphs. While the simulations testing the effects of 
the explosion energy and the binary separation were carried out 
with 2x10^ particles in the companion star, the study based on 
the realistic progenitor system structure was set up with 2x10^ 
and 4x10^ particles representing the companion. Therewith, the 
mass of the single particles and, as all particles had the same 
mass by construction, the number of particles representing the 
supernova explosion was fixed. 
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Fig. 3. Stripped mass for different supernova energies in model 
rp3^0a. 



The stripped mass as a function of the supernova energy is 
shown in Fig.|3] The relation is linear in good approximation and 
can be fitted by 



^stripped = 2.6 X 10 



-1 -^kir 



,SN 



lO^'erg 



M„, 



(2) 



assuming that the offset is zero (without this constraint the offset 
is only 1.3 x 10"^). We emphasize that, although the functional 
form of the relation may be generic, the particular parameters 
of the fit apply to model rp3_20a only. For different companion 
structures, for instance, the values of the parameters are expected 
to change. This should be kept in mind for the fits presented 
below as well. 

Note that the studied energy range covers only a factor of 2 
and therefore the stripped mass also changes only by a factor of 
2. The energy of the W7 model of 1.23 x 10^' erg corresponds to 
an intermediate case. These results indicate that changes in the 
supernova energy have only a small effect on the stripped mass 
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Table 2. Paramete rs of the progenitor mod els and of the companion star after interaction with the supernova ejecta. 

Values taken from ilvanova &■ Taarm2004) : masses of the companion star at the beginning of the mass transfer Md,i and at the time 

of the explosion Md,f , length of the mass transfer period Aftr, orbital period Pf, distance between white dwarf and its companion star 

just before the explosion Of ; 

Results of the simulations: mass stripped from the companion Mstripped and its kick velocity yyck, 5000 s after the explosion. 



as compared to the other parameters discussed below. The kick 
velocity can be fitted by a power law as 



yjdck = 2.5 X 10 
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erg 
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A simple argument for this behavior may be that velocity of the 
companion is mainly given by the momentum exchange from 
the supernova ejecta that hit the star As the velocity of the ejecta 
increases with the square root of the kinetic energy, the velocity 
of the star should increase accordingly. The exponent in eq. (O is 
only slightly larger than 1 /2 and thus this simple picture captures 
the process rather well. 

4.2. Distance 

In order to test the influence of the separation, we again choose 
model rp3_20a (see Table |2] for the parameters). Here, the dis- 
tance d between WD and companion was varied in the range 
of 1.5 ... 3 X 10" cm. All other parameters were kept constant. 
The relation between distance and stripped mass for this model 
is shown in Fig. |4] From the largest to the smallest distance the 
stripped mass increases by a factor of 10. The relation follows a 
power law in good approximation and can be fitted to 



M,s,ripped = 2.38 X 10^^ I — 
\cm 



-3.49 



Mo. 



(4) 



The kick velocity depending on the separation d also follows a 
power law. It is approximately given by 



wkick = 2.5 X 10 



23 



cms 



(5) 



The fraction of the supernova ejecta that hits the star scales 
with the inverse square of the distance treated as an independent 
parameter as discussed above. Remembering that the stripped 
mass seems to scale linearly with the explosion energy, one may 
expect a geometrical scaling as Mstripped <^ d^'- This does clearly 
not fit our results. It indicates that the connection is more com- 
plex, and the linear correlation between explosion energy and 
stripped mass is possibly only a result of the small energy range 
tested. In addition the simple scaling neglects that the companion 
star is not flat disc but a sphere. The whole process also seems to 
depend on the time evolution of the density of the ejecta hitting 
the star. 
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Fig. 4. Stripped mass versus binary separation for model 
rp3_20a. 



4.3. Companion star 

Table |2] shows the stripped masses and the kick velocities of 
the remnants for the six dif ferent progenitor systems suggested 
by llvanova & TaamI (|2004|) . The supernova model used in all 
models is the original W7 model with a total kinetic energy of 
1.23 X 10^' erg. The variations in the stripped mass between the 
2x10^ particles and the 4x10'' particles simulations were less 
than a few percent indicating numerically converged results. 

The stripped masses range from 0.01 Mq to 0.06 Mq for 
the different setups. This is si gnificantly less than previous re- 
sults r eported by other authors dlVIarietta et al.ll2000t llVIeng et al.l 
l2007h . Compared to the 0.15 Mq result of llVIarietta et al.1 (l2000h 
we find a factor of 3-15 less stripped material. 

This deviation is attributed to the binary evolution of the pro- 
genitor we take into account in the present study. The binary evo- 
lution significantly aff'ects the properties of both the companion 
star and the geometry of the binary configuration. 

The main effect on the companion star is illustrated in Fig.|5] 
where the density profile of the companion star in model rp3_20a 
at the time of the explosion is plotted in comparison with the 
density profile of a star in thermal equilibrium that was evolved 
as a single object to the same mass and nuclear age. The equi- 
librium star features a much larger radius than the binary com- 
panion. This is a result of the mass transfer phase in the binary 
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system, that removes the outer layers and leads to a more com- 
pact star. The mass loss occurs rather fast: the duration of the 
mass transfer phase is about a factor of 10 less than the Kelvin- 
Helmholtz time-scale of the stars in our models. Therefore, the 
star is not able to adjust its structure to the loss of the outer layers 
and it shrinks into a more compact object. Finding less stripped 
mass when taking into account the modification of the compan- 
ion structure due to the binary evolution is thus not surprising. 

In order to corroborate this interpretation, we set up model 
rp3_20a with a companion structure that corresponds to that 
of a star of equal mass in thermal equilibrium. In this model, 
0.066 Mq are stripped from the companion by the supernova. 
This is a factor of 2 more than for the original rp3_20a model. 
The diff'erence can be explained by the larger, less bound en- 
velope of the equilibrium star that can be stripped away more 
easily. Moreover, the larger radius causes an extended interac- 
tion area that also leads to a slightly larger kick velocity of 
51.2kms ' for the equilibrium star companion model evolved 
in isolation as compared to the 46.6 km s"' of the binary star 
companion model. 

Additionally, in the model we use the companion star is 
slightly more massive (1.17 M© at t he time of the explosio n) than 
the solar mass companion star of [Marietta et al.l (l2000l) . This 
leads to a slightly larger radius of the equilibrium star compan- 
ion and therewith a larger separation distance and a decreased 
stripped mass. 

These effects together explain the sig nificant difference be - 
tween our results and the previous work of lMarietta et al.l (l200d) . 
who used a companion with the structure of an isolated equi- 
librium star We note, howev er, that our interpre tation is not in 
agreement with the work of iMeng et alj ( l2007h . who found a 
lower limit on the stripped mass of 0.035 Mq taking into account 
binary evolution. This is probably be due to oversimplifying as- 
sumptions they made in their analytical treatment of the ejecta- 
companion star interaction. 

The kick velocities of the companion star after the impact of 
the supernova ejecta we find in our models vary from 17 km s ' 
to 61 kms"^ (see Table|2]). The velocities roughly increase with 
the size of the companion stars (as the cross section increases) 
and decrease with larger separation distances (as the transferred 
momentum decreases) and masses (for the same transferred mo- 
mentum) of the companion stars. 
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Fig. 5. Comparison of density profiles of the companion star in 
model rp3^0a at the time of the explosion and a single star with 
the same mass and nuclear age. 



5. Observational implications 

Observations that may help to constrain the nature of the SN la 
progenitor systems are the detection of hydrogen in the spectra 
of these events, signatures in the spectrapolarimetry data due to 
an asymmetric morphology of the ejecta as a result of the inter- 
action with the companion star, and a direct observation of com- 
panions in supernova remnants. We discuss these possibilities in 
the light of the presented study. 



5. 1 . Hydrogen detection in SN la spectra 

Our results are consi stent with the c onstr aints on hydrogen i n 
the ejecta as given by iLeonardI ( l2007h and lMattilaet al.l ( l2005h . 
However, none of our models has strip ped hydrogen mass far be- 
low the upper limit of iLeonard ilOOJh of 0.01 Mq. As these Um- 
its result from a non-detection of hydrogen in nebular spectra, 
our simulations predict hydrogen detection not far below these 
values. Conclusions on the validity of the single-degenerate sce- 
nario depend on whether the systems studied here are represen- 
tative for Chandrasekhar mass WD+MS SN la progenitors and 
whether the observed events fall into this class. If both were true, 
lowering the observational upper limits of hydrogen in the ejecta 
of SNe la by another order of magnitude would exclude this pro- 
genitor scenario. 

Yet it i s important to n ote t hat a quite sirnple mo del was em- 
ployed bv lLeonardl(l2007h and iMattila et al.l J2005h to constrain 
limits on the hydrogen mass from the observations. A more rig- 
orous approach would be to use the results of hydrodynamical 
simulations such as presented here as an input for full radiative 
transport calculations. From these calculations it will be possible 
to predict whether hydrogen lines should be visible in the spec- 
tra and how strong they should be at a given epoch. This issue 
will be addressed in a forthcoming study. 



5.2. Hole in the ejecta 

The supernova ejecta do not only affect the companion star, but 
are also affected themselves by the impact. Fig.|6]shows the ma- 
terial after the impact. The left column shows both the supernova 
ejecta and the companion star material including the hydrogen 
stripped from the star The right column shows only the mate- 
rial that was part of the supernova ejecta at the beginning of the 
simulation. 

The supernova ejecta that were spherically symmetric in the 
beginning are clearly asymmetric after the impact. In the wake 
of the companion star, a cone-like hole in the supernova ejecta 
is visible. To some degree it is filled with material that is ejected 
from the companion star. At the borders of this hole, the super- 
nova ejecta are slightly denser, because the material missing in 
the hole was transfered there. The opening angle of the cone- 
like hole is about 45° (see the top row of Fig. |6] showing the 
ejecta 1000 s after the explosion). T his result is consistent with 
the findings of lMarietta et al.l ( |2000|) . However, the area in which 
the supernova ejecta are affected by the impact on the compan- 
ion star is as large as 90° (see middle and the lower row of Fig. |6] 
showing the density of supernova ejecta after 3000 s and 5000 s). 

A more detailed view of how the supernova ejecta are mixed 
with the material stripped from the companion star is given in 
figure |2| It shows for the the same setup the relative amount of 
material that originally belonged to the companion star with re- 
spect to the total amount of material. 

Qualitative ly, the structure of t he ejecta resembles the simple 
model used bv lKasen et al.l (|2004|) to explore the effect of a hole 
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Fig. 6. Snapshots of the evolution of the supernova ejecta in the rp3_20a scenario. The left column shows supernova ejecta and the 
companion star In the right column all material of the companion is cut out, leaving only the supernova material. The plots use 
cylindrical coordinates. The radial coordinate is averaged over angle. Color-coded is the density. 



in the ejecta on spectra and luminosity of SN la. With it he was 5.3. Identifying companion stars in supernova remnants 

able to reproduce observed spectrapolarimetry observations of 

cj^ T The remaining companion star should have a different velocity 

than its surrounding stars. This velocity is determined by the or- 
bital velocity of the star at the moment the white dwarf explodes 
and of the kick it gets from the impact of the supernova ejecta. 
The former is perpendicular to the connecting line between the 
white dwarf and the companion star. For our models it ranges be- 
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coordinates with the radial coordinate averaged over angle are used. 



tween 130 kms^' (rp3^0c) and BSOkms^' (rp3_28a). The latter 
velocity is a result of the impact of the ejecta and is therefore 
aligned with the connecting line. In our simulations it reaches 
values from 17km s"' to 50 km s"' . Therefore, the orbital veloc- 
ity clearly dominates the velocity of the companion star relative 
to the center of the supernova remnants. The velocity of the star 
Tycho G whi ch was identified as prog enitor of Tycho Brahe's 
supernova bv lRuiz-Lapuente et alj ( l2004i) features a modulus of 
its spatial velocity of 136kms '. It thus falls into lowest part of 
the range we find in our simulations. The properties of Tycho G 
are thus consistent with the predictions of our models. 



6. Conclusions 

We studied the impact of the ejecta of SNe la on main se- 
quence companion stars in the context of the single-degenerate 
Chandrasekhar-mass scenario with hydrodynamical simulations. 
To this end, the (cosmological) GADGET2 smoothed particle 
hydrodynamics code was adapted to the stellar problem and em- 
ployed in numerical simulations. It was shown that this SPH- 
based approach is capable of reprodu cing previous results ob- 
tained with a grid-based 2D scheme bv lMarietta et al.l (l2000l) . A 
resolution study indicated that with a few million particles the 
simulations yield numerically converged results. 

We showed that the mass stripped from the companion star 
by the impact of the ejecta depends on the their kinetic energy 
and the binary separation. While the latter affects the mass of the 



stripped material significantly, the former has only a minor effect 
on it. This is due to the fact that the supernova explosion energy 
can only vary in a relatively narrow range given the restricted 
amount of fuel available for the nuclear energy generation. 

The SPH approach was used to analyze the impact in a 
number of more realistic progenitor models than those em- 
ployed in previous studies. For these, the companion stars 
were constructed with the stellar evolution code GARSTEC 
mimicking binar y mass transfers with the parameters given by 
llvanova & TaamI (12004!). In the hydrodynamical impact simu- 
lations, we found about one order of magnitude less hydrogen 
material stripped off the companion by the impact of the super- 
nova than predicted by previous studies. The main reason for 
this difference is a modified, more compact stellar structure of 
the companion star in combination with a resulting variation in 
the separation distance of the progenitor system. In particular the 
more compact state of the companion impedes the mass loss in 
the impact. 

The reduced amount of hydrogen mixed into the ejecta 
of the supernova as predicted by our simulations leads to an 
agreement with observational studi es of SN la nebular spectra 
JMattila et al.1 120051: [Leonard! l2007l) . This removes the former 
disagreement between the available observations and simula- 
tions of the WD+MS progenitor system. Thus, to current knowl- 
edge, such a progenitor scenario is admissible in the context dis- 
cussed here. 
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However, since the hydrogen masses predicted by our simu- 
lations are not far below the current observational upper limits, 
it may be possible in the near future to confirm or reject the stud- 
ied progenitor scenario by either detecting hydrogen in SNe la 
or lowering the limits by another order of magnitude. A stringent 
way of analysis would be to calculate synthetic spectra directly 
from the presented simulations and to compare the results with 
observations. This will be tackled in a forthcoming study. 
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